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Abstract: An approximate solution of the Klein-Gordon equation for the general Hulthen-type 
potentials in D-dimensions within the framework of an approximation to the centrifugal term is 
obtained. The bound state energy eigenvalues and the normalized eigenfunctions are obtained in 
terms of hypergeometric polynomials. 
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O , I. INTRODUCTION 

^\ 

\^ ' The search for exact solutions of wave equations, whether non-relativistic or relativistic, has been an important 
] research area since the birth of quantum mechanics. The generalized Hulthen potential [l|-[18| is given by 

V{r) = -Za (1) 

■ 

I where a is the screening parameter and Z is a constant which is identified with the atomic number when the potential 
. is used for atomic phenomena. The Hulthen potential is one of the important short-range potentials which behaves 

I ^ I ' like a Coulomb potential for small values of r and decreases exponentially for large values of r. The Hulthen potential 
has received extensive study in both relativistic and non-relativistic quantum mechanics [ll|-[l8l|. There is a wealth 

O^l ' of literature on the use of the Hulthen potential as an approximation of the interaction potential in a number of 
[ area in physics such as nuclear and particle physics [H, atomic physics [ll-Q, solid state physics [H] and chemical 

■ physics [6|, see also and the references therein. uSortunately, quantum mechanical equations with the Hulthen 
potential can be solved analytically only for states with zero-angular momentum [H-[ll]. Recently some interesting 
research papers [iT'l-fls] have appeared to study the /-state solutions of the relativistic Klein-Gordon equation with 
Hulthen-type potentials. The main idea of their investigation relies on using the approximation of the centrifugal 
term p- by means of p- « (^■^_^-ary2 ■ Their results show that this approximation is in good agreement with the 

t — ■ other methods for small a values. The purpose of the present work is twofold: (i) to extend the Z-state approximate 
solutions of Klein-Gordon equation with the generalized Hulthen potentials to arbitrary dimension; (ii) to 

^ compute the normalization constant of the approximate wave functions that seems to have been overlooked by many 
researchers fH-flsj. 



II. THE KLEIN-GORDON EQUATION IN D DIMENSIONS 

The £)-dimensional Klein-Gordon equation for a particle of mass M with radially symmetric Lorentz vector and 
Lorentz scalar potentials, V{r) and S{r), r = ||r||, is given (in atomic units h — c=l) [l^[20| by 

{-Ai5 + [M + 5(r)]2}*(r) ^ [E ~ V{r)]H{v), (2) 

where E denotes the energy and A/j is the /^-dimensional Laplacian. Transforming to the D dimensional spherical 
coordinates (r, 9i . . . Ojj^i), the variables can be separated using 

'fir) ^R{r)Yi,_,,..jM-- -00-1) (3) 
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where R(r) is a radial function, and Yi^^ -^^,,, ^1-^(61 . . . 0_d_i) is a normalized hyper-spherical harmonic with eigenvalue 
l{l + D — 2), I — 0, 1, 2, ... . Thus, we obtain the radial equation of Klein-Gordon equation in D dimensions by 
substituting Eq.® into Eq.^ 



(f D-1 d l{l + D-2) 
dr"^ r dr 



+ [M + S{r)Y ~[E~ V{r)Y \R{v) = 0. (4) 



Writing R as R{r) = r ^)/^w(r) gives 

- ^ + ' ' + [W + S{r)r -[E- V{r)r] \u{r) = 0. (5) 



where k = D + 21 and u{r) is the reduced radial wave function satisfying w(0) — 0. 



III. BOUND STATES FOR GENERALIZED HULTHEN POTENTIAL 

We consider the vector and scalar Hulthen potential defined as, 

nr)=-^^^^, S{r) = -^^^, (6) 
1 — qe 1 — qe "'^ 

where Vq and 5*0 are the depth of the vector and scalar Hulthen potential respectively and a is the screening parameter 
and g 7^ is the deformation parameter. Substituting ([6]) in the Klein Gordon equation (O, we obtain 

, (fc-l)(fc-3) i2MSo + 2EVo)e-"^ (Si-V,^)e-^"n ^ ^ , . 



dr'^ 4^2 I _ q^-ar (^i _ qf^-ary 



u{r) = [E^ - Ar)u{r). (7) 



In order to obtain analytic solutions of this equation, we have to use an approximation [l7j-p^ for the centrifugal 
term similar to that used for the non-relativistic cases. We, thus, follow IlH-llll and use k. , "^'^^"^^ for 

r \L qe ) 

the centrifugal term. This approximation is valid for q = 1, however, we follow the model used by Qiang et al [iSj . 
This allow us to write Eq.© as 



dr2 1 - ge-"'^ (1 - qe-"'')^ 

Eq.® can be further simplified using a new variable z — ge""'' (r e [0, 00), z E [q, 0)) 

/?! i(/c-l)(fc-3) f32 



dr'^ z dz z^ z(l — z) qz{l — z)^ (1 — z)^ 

where we used the dimensionless parameters given by 



u{z) = 0, (9) 



VM^ - E^ 2iMSo + EVo) - Vp' 

e= , pi — 5 , P2 = o n ■ UUj 

a a'^q a'^q^ 

We now look for a solution of ^ in the form 

u{z)^z^{l-zYf{z). (11) 

In this case, Eq.® reads 



iPi - (2e + l)S-5^+S-P2)il-z) + 5^-S-fh-i^{k-l){k- 3) 



z(l-z)2 



fiz) = 

(12) 
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We may now choose 6 such that 

1 



which yields 



6^-S-(32-^{k-l){k-3) = (13) 



1 1 



5 = S± = -± ^ ^^'(l + m + q{k - l)(fc - 3) (14) 



where (5 = (5+ for g > and 5 = &^ for g < 0. Thus Eq. p^ reduce to 

This equation is a special case of a more general differential equation discussed in (29j , namely 



2az^+i 2(m + l)\ 



AT 



r(^) - (^ ^ „ ^^^,+2 -^-^ jn^) - _ ^ = -1,0,1,... (16) 

which has exact solutions, for n = 0, 1, 2, . . . , given by 
if 

u;„(iV) = n(Ar + 2){{n{N + 2) + (2m + l))b + 2a). (18) 

Here C„ is the normalization constant and 2F1 (a, 6; c; n) is a special case [l^l of the generalized hypergeometric 
function 

re ^ (Q^l)n • • • («p)n ^" /-ifiN 
pF,(ai, ...,ap;ci, ...,c,;z) = 2^ — (19) 

\Cl)n ■ • ■ [C-qjn n. 

where the Pochhammer symbol (a)„ is defined by (a)„ = T(a + n)/T{a). By putting = —1 and 6 = 1 in pop . using 
P^ . with a — S and w„(— 1) = n(n + 2e + 2(5), we obtain the energy spectrum of © as 



- - (1 + 2n)5] ~ \{k - l){k - 3) 
2q{n + d) 

_ q[Pi -in + 6)^] + q[S^ - ^ - ^ (fc - 1) (fc - 3)] 
2<7(n + 6) 

= 2(;rM)^2(" + ^) 
where we have used (|13p . Furthermore, the exact solutions of ([72]) . using (flTl) . now reads 

/„(z) = (-l)" C„ (l + 2eW)„ 2-Fi(-n,26W+2,5 + Ji;l + 2eW^^) (21) 
Therefore, we can write the total radial wave function (fTTj) as follows 



,{z) = C„ z^"'(l-z)Si^i(-n,2el^-)+25 + n;l + 2eW,;z) 

= C„z^S^^'(l-z)^/^l^^"'^^^--^)(l-2z) (22) 



where we used the definition of Jacobi polynomials [3l| given by 



(23) 
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IV. NORMALIZATION CONSTANT 



In this section, we compute the normalization constant C„ appear in (|22p . As far as we aware, the normahzed 
energy eigenfunctions have not been exphcitly worked out in the hterature for this case. It is straightforward to note 
that the normahzation constant C„ can be computed using |i?^(r)|r'^^~^^(ir — \u{r)\'^dr — — J° = 1 

because of the substitution z — qe^"^ . Thus, we have by mean of (I22p that the normalization constant C„ is given by 



Ik) - r ,, . ,, , 1 2 



/ z^^" -1(1-^)2* 2i^i(-n,2eW +2,5 + n;l + 2eW dz ^ a. (24) 







Using the series representation (jl9p of the hypergeometric function 2^^!, being a polynomial of degree n in z, we have 



i=0 j=0 







The definite integral in (|25|) is just the integral representation of Incomplete Beta function [32| . 

Bq{x,y)^ r e-^{l~t)y-^dt, 3fJ(a;) > 0, 3fi(y) > 0. (26) 
Jo 

Therefore, Eg. ([^5]) now reads 

(-^M2^^+ 2. + n)a-n),(2^) + 2. + n), + , + 2, + l) = « (27) 

t^,t^ (l + 2.W).^! (l + 2.«),j! 

On integrating by parts of (|26p . we can find that the Incomplete Beta function satisfies the recurrence relation 

B,{x + l,y) = -^B,{x,y)-'^^^^^-^. (28) 

X + y X + y 

which in turn can be written in terms of the normalized version of the Incomplete Beta function Iq{x,y) — 
Bq{x,y)/B{x,y) as 



= Iqi^- 






g^-i(l - 


q)y 








[x 


~ l)B{x 










= Iq{x- 


2,2/)- 




g^-2(i _ 


q)y 


q^-i(l - 


q)y 




(x 


- 2)B{x 


-2,2;) 


{x-l)B{x 


-i,y) 




= Iqi^- 


3,2/)- 




g^-3(l - 


q)y 


q^-2(l - 


q)y 


9=^-1(1 - q)y 


{x 


- 3)B{x 


-3,y) 


{x-2)B{x 


-2,y) 


{x - l)B{x - l,y) 



/,(.-^,,)-,-(l-,)^g ^^_^^^^^_^^^^ , ^ = 1,2,... (29) 



Thus 



/,(2eW+^ + J,2^ + l) = /,(2.W, 25 + 1) - g2^^^'+^+^(l - 5)^^+1 -^^^ — ^I^^ 

tr[ {2er k)B{2er +i+j-k,2S + l) 

(30) 

which allow us to compute the Incomplete Beta function in ^27]i as Bq{2eh^ + i + j, 2(5 + 1) = 

B{2e^n'' + i + 3,25 + l)/^(2el''^ + i + j,26 + 1) for i,j = 0, 1, 2, ... . Note, in the case of i = j = 0, the sum 
in (f30|) is equal to zero. 
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Although the discussion above for computing the normahzation constant assumed that q G (0, 1), the computation 
for arbitrary q can be performed using the analytic expressions 32]: 

q^il-q)y-^^{l-y)kf q 



^^^i^^.^r(M-.;i + .;^) (31) 



for q e (-00, 0) U (0, i), and 



R( ^ R( ^ g--^(l-g)^ ^ (l-x)fc ^g-l^' 



= B{x, y) - ^l^ii^^Fi (1, 1 _ x; 1 + (32) 

y q 

for (? > 5- In the case of g = 1, i.e. the Hulthcn potential, Eg. ([5^ yields Bi{x,y) = B{x,y), consequently, Eg. ([^5]) 
becomes 

y^y. (-^).^(2^^+ 2. + ^)^(-»M2^^ + 2. + .),^ + . + 2. + 1) = « (33) 
Uf^o (l + 2^^).^! (1 + 2.W),,! 

Using the definition of Beta function 's?] in terms of Gamma function B{x,y) — ^Y(^^^y write ([33| as 

^2 (-n)»(2£i^^ +26 + n), i-n)j{2e'n^ +26 + n), r(2£l^^ + i + j)V{25 + 1) ^ ^ 

By means of the definition of Pochhammer symbols (a)„ = '"p"^"-* , we have 

^2 r(2£i''^)r(2,5+l) i-nU2e'r!'^ + 26 + nU2€i% (-n)j(2ei''^ + 2^ + n)j(2ei''^ + z)^ 
" r(2el"-) + 2<5 + 1) (1 + 2e^^^),(2e(^) + 2<5 + 1),*! ^ (1 + 2ei''>),{2eit^ + i + 26 + 1),.?! 

Thus, by using the series representation of the hypergeometric series 3F2, again Eg. p^ . Eg. ([55)) then reduce to 

^2 ^^ hn),{2ei'H26 + nU2ei^'>), ( -n, 2el^^ +26 + n, 2ei^^ +i \ _ a 



which can be used to compute the normalization constant for n = 0,1,2,.... In particular, for the ground-state 
n = 0, we have 



Co^J ^ ^ -. (37) 

V S(2e[,'=\25+1) 



V. CONCLUSION 



In this work, we have extended the approximate analytic solutions of Klein-Gordon equation with vector and scalar 
generalized Hulthen potential to arbitrary dimension D. The analytical energy equation and the normalized radial 
wave functions expressed in terms of hypergeometric polynomials are given. When D = 3, our results normalize the 
approximate analytic solution for bound states obtained in [T^ and for Klein-Gordon equation with generalized 
Hulthen potentials for nonzero angular momentum. 
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